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Abstract 

We show that sparse spherical harmonic expansions can be efficiently recovered from a small 
number of randomly chosen samples on the sphere. To establish the main result, we verify the 
restricted isometry property of an associated preconditioned random measurement matrix using 
recent estimates on the uniform growth of Jacobi polynomials. 

1 Introduction 

Compressed sensing has triggered significant research activity in recent years. It predicts that sparse 
signals can be recovered from what was previously believed to be highly incomplete information. 
In this work, we show that functions on the sphere S 2 = {x G M 3 , ||x||2 = 1} that have a sparse 
or compressible representation in the spherical harmonic basis can be recovered from a number 
of samples that scales (essentially) linearly with the sparsity level. This can be viewed as an 
extension of existing results |B] for sparse recovery of trigonometric polynomials on the circle. 
Since the (L2-normalized) spherical harmonic basis functions are not uniformly bounded, standard 
compressed sensing theory does not apply directly to sparse recovery on the sphere. Instead, 
we appeal to recent results in [1] concerning sparse recovery in orthonormal polynomial systems. 
Since orthonormal polynomials blow up sufficiently quickly and uniformly at the endpoints of their 
domain, one preconditions in order to transform the problem into that of sparse recovery in a 
uniformly bounded system. The decomposition of spherical harmonic basis functions into tensor 
products of trigonometric polynomials and Jacobi polynomials allows to prove the main result: 
any degree- D polynomial on the sphere (that is, with N = D 2 coefficients) consisting of at most s 
spherical harmonic basis elements can be efficiently recovered from m ~ 

sATi/4 j og 4 independent 

sampling points drawn uniformly with respect to a certain measure (see below). We establish this 
result by verifying the restricted isometry property (RIP) of an associated random matrix. 

2 Background and notation 

The spherical harmonics Yg, —l< k < £,k > form an orthonormal basis for the Hilbert space 
of square-integrable functions on the sphere. They are orthogonal with respect to the spherical 
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surface measure Q. In spherical coordinates (4>,6) G B := [0, 7r] x [0, 27r), (x = cos(0) sin(^>), y = 
sin(0) sin(^), z = cos(c/>)) G S 2 , this orthogonality relation becomes 

/ 27r r e)Yf($J) S m((f>)d^de = Su>5 w . (i) 

io io 

Here, z denotes the complex conjugate and 5 mn is the Kronecker delta. The spherical harmonics 
may be expressed as 

y/(0^) = e^(sin0)IS[ fc i| fc |(cos^), ((f), 6) G B, (2) 

where the Jacobi polynomials (p")^L with parameter a > — 1 are the orthonormal polynomial 
basis on the interval [—1,1] with respect to the measure dv(x) = (1 — x 2 ) a dx In particular, 
the Lebesgue measure a = generates the Legendre polynomials, while the Chebyshev measure, 
a = —1/2, generates the Chebyshev polynomials. 

Spherical harmonic expansions on the sphere are analogous to Fourier series expansions on the 
circle. Functions on the sphere of the form 

D-l i 

1=0 k=-t 

are called harmonic polynomials of degree D — 1. Note that N = D 2 spherical harmonic basis 
elements generate harmonic polynomials of degree D — 1. We will call a harmonic polynomial s- 
sparse if its coefficient vector c = (q^) G C n has cardinality at most s; i.e. ||c||o := |{(^, k) : C£ t k / 
0}| < s. More generally, the degree to which a harmonic polynomial can be well-approximated by 
its s most significant coefficients can be quantified using the concept of best s-term approximation 
error which is defined, for a vector z G by 

o- s (z)i= inf [|y-z||i. 
v-\\y\\o<s 

We say that a harmonic polynomial fcty is compressible if <r s (c)i decays quickly as s increases. 



3 Main results 



We aim to recover sparse harmonic polynomials on the sphere from only a few function samples. 
Note that samples yj = g(4>j,6j), j = 1, . . . , m, of a D — 1-degree harmonic polynomial g(4>, 9) = 
lLje=o Ylk=-£ c £,kY£((j), 0) may be expressed concisely in terms of the coefficient vector c = (c^j.) G 
according to 

y = $c. (4) 
where $ is the m x N matrix defined component-wise by 

$ Uk , i) =Y £ k (ct> j ,e j ). (5) 

We are interested in solving the system of linear equations Q in the underdetermined setting 
m < N, and in particular, to single out the original sparse coefficient vector c from among the 
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infinitely-many solutions. The compressed sensing literature has suggested various reconstruction 
algorithms for sparse recovery; for simplicity we focus only on ^-minimization [7\ in this paper. 

The Loo-norm of the spherical harmonics Y^{4>,6) increases with the degree £ according to 
sup_ e<k< i ||Y^||oo = sup ||l^°||oo = £, and this extremum is obtained at the spherical caps (j) = 0, ir. 
This means that the linear system Q is in general ill-conditioned. It is known (see Proposition [6]) 
that | (sin 0) 1 1 ' 2 'Y^(<f>, 9)\ ~ consequently, we precondition the system Q for numerical stability, 
multiplying both sides by the m x m diagonal matrix A with entries Aaa = (sin^) 1 / 2 , 

Ay = A$c. (6) 

Our main result is that any s-sparse harmonic polynomial on the sphere of maximal degree 
D can be recovered efficiently from a number of samples m that scales linearly with the sparsity 
level and sublinearly with the degree. This reconstruction is moreover robust with respect to noisy 
samples and passing from sparse to compressible vectors. 

Theorem 1. Let m,s, and N be given integers satisfying 

m>Cs log 3 (s)iV 1/4 log(JV). (7) 

Suppose that m coordinates on the sphere (</>i, 6i), . . . , ((fim, 9 m ) are drawn independently from the 
uniform measure on B = [0, 7r] x [0, 2tt) . 

Let be the m x N spherical harmonic matrix ^ and let A& be its preconditioned version 
With probability exceeding 1 — _/V~ 7log ^ the following holds for all harmonic polynomials 
g(4>,9) = J2e > =o 1 J2k=-e c i,k'Yf k {4',0). Suppose that noisy sample values yj = g(cf)j,9j) + r/j are 
observed, and that ||?7||oo 

< e. Let 

c = argmin ||z||i subject to \\A$z — Ay\\2 < y/rne. (8) 

Then 

n , Cia s (c)i 

V s 

The constants C,C\,C2, and 7 are universal. 

We expect that the bound ([7]) is not optimal, but so far we have not been able to remove the 
polynomial factor N l l A . 

4 Numerical experiments 

In Figure 1, we plot two phase diagrams illustrating the success of ^i-minimization (i.e., Equation 
^ with e = 0) in recovering s-sparse harmonic polynomials g(<ft,0) = X^^o 1 Ylk=-e c e,kY^((f), 6) 
from m sample values g(4>j, 0j). In each plot, we fix the degree to be D = N 1 / 2 = 16 and vary the 
sparsity level s and number of measurements m. We form an s-sparse coefficient vector c = (c^fc) 
by choosing a random support set from Nfl [l,iV] of cardinality s and assigning independent and 
identically distributed Gaussian weights as the coefficients to this support; we then draw m sampling 
points which we use to recover the s-sparse vector using £i-minimization. For each pair (s,m), we 
record the frequency of success of £i-minimization out of 20 trials. In Figure 1(a) the sampling 
points are chosen independently from the product measure d(j>d6 on B, which has higher sampling 
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Figure 1: Phase transition for sparse recovery on the sphere 
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density around the spherical poles and for which the sparse recovery results of Theorem 1 apply. We 
observe a large region of phase space (in black) corresponding to uniform recovery, similar to that 
observed in the phase transition curves obtained for other compressive sensing matrices |5j . On the 
other hand, in Figure 1(b) the sampling points are chosen independently from the uniform surface 
measure sin (<fr)d(f)d6. In this case, ^i-minimization fails to recover sparse harmonic polynomials for 
essentially all parameters (s,m). 



5 Sparse recovery via Restricted Isometry Constants 

We prove Theorem 1 by showing that the preconditioned spherical harmonic matrix <I> = in 
^ satisfies the restricted isometry property (RIP) [7]. 

Definition 2 (Restricted isometry constants). Let ^ G C mxN . For s < N, the restricted isometry 
constant 5 S associated to \I/ is the smallest number such that 

(l-<5.)l|c||2<ll*c||l<(l + * a )||c||l (9) 

for all s-sparse vectors c G C . 

Informally, we say that the matrix ^ "has the restricted isometry property" if 5 S is small for 
s reasonably large compared to m. For matrices satisfying the restricted isometry property, the 
following ^-recovery results can be shown [3 [8]. 

Theorem 3 (Sparse recovery for RIP-matrices). Let * G C mxN . A ssume that the restricted 
isometry constant of ^ G ^mxN sa f{ s -p es 

5 2s < 3/(4 + V6) ~ 0.4652. (10) 

Let x G C N and assume noisy measurements y = + rj are given with || 77 1| 2 < s. Let x& be the 
minimizer of 

arg min \\z\\i subject to \\$>z — wlU < e. (11) 

zeC N 
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Then 

\\ x - x * h < Cl ^^ + C,:. (12) 



for some constants C\,Ci > i/iaf depend only on &2s- In particular, if x is s-sparse then recon- 
struction is exact, = x. 

A general setup for matrices having the restricted isometry property are those associated to 
bounded orthonormal systems [21 E] - 

Theorem 4 (RIP for bounded orthonormal systems). Consider an orthonormal system of func- 
tions ipj, j £ [1,N] n N on a measurable space Ai endowed with a probability measure v, that is 
Im ^j^kd-v = Sjfi. Consider the matrix ^ G £ r mxN w itfa entries 

®l,k = M x t), i € [1, m] n N, k e [1, N] n N, 

formed by i.i.d. samples X£ drawn from the measure v. Suppose this system has the uniform bound 
K = sup je[Ar] Halloo = sup ie[Ar] sup xeC \tpj{x)\. If 

m > C5~ 2 K 2 s log 3 (s) log(iV), (13) 

then with probability at least 1 — jV -710 ^ 5 ), the restricted isometry constant S s of satisfies 
5 S < 5. The constants C > and 7 > are universal. 

An important special case is the matrix associated to samples of the trigonometric system 
(exp(27rijx))^ ) 1 chosen from the uniform measure on [0, 1], which has the optimal uniform bound 
K = 1. Another example is the sampling matrix associated to the Chebyshev polynomial system. 
In this case, K = s/2. 

6 Sparse recovery in spherical harmonic systems 

Recall from ([2]) that the spherical harmonics can be expressed as tensor products of complex 
exponentials in 9 and orthogonal polynomials in cos(</>). Since the latter are not uniformly bounded, 
spherical harmonics do not fall directly into the scope of bounded orthonormal systems. To get 
around this obstacle, we proceed in a similar fashion to p], and use estimates on the uniform rate 
of growth of orthogonal polynomials in order to precondition the spherical harmonic system. 
First we will need the following growth estimates. 

Proposition 5. Consider the weight function v(x) = (1 — x 2 ) a on [—1,1], and let (p%) n be the 
associated orthonormal polynomial system. Then, for all x & [—1,1], the following holds. 

1. If a = 0, the associated polynomials (p^) are the Legendre polynomials and satisfy 

{l-x 2 )^\pl(x)\<2/^. (14) 

2. For any a > 0, 

(l-X 2 )V 4 +«/2| p a( x )| < Ca ( 15 ) 
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3. Ifa> 3, then 

(l-x 2 ) 1/4+a/2 |KWI<Ca^( 1 + ^) 1/12 , (16) 
where C > is a universal constant. 



The bound (14) for Legendre polynomials is classical and known to be tight, and the more 



general bound (15) is also classical; see [1] for more details. The more refined bound (16) was 
derived only recently in [3]. Although the result is stated in [3] for the parameter range a > 1, n > 6, 
it is not hard to verify that the key estimate (Lemma 8 in [9]) is valid also for the parameter range 
a > 3,n > 0. 

Using these bounds in conjunction with the tensor product representation Q, we arrive at the 
following rate of growth for the spherical harmonics. 

Proposition 6. For all £ 6 Z + and — £ < k < £, 

(sin 4>) l l 2 \Yf(<j), 0)| <C(^ + 1) 1/4 VxG [-1,1]. 
We can now state the proof of Theorem 1. 
Proof of Theorem 1. 

Consider, for < £ < D — 1, —£< k < £, the functions 

Q k e (<P,0) = ( S m^ 1 / 2 Y e k ( ( f ) ,e). (17) 

By Proposition [6j 

sup HQ^Iloo < CiV 1 / 8 

0<i<VN-l-£<k<£ 

for a universal constant C. Because the spherical harmonics Y k are orthonormal with respect to 
the uniform measure shx{(j))d(j)d9 , the Q^'s are orthonormal with respect to the product measure 

2 " f * Q k t(<t>,6)Q%((l>,6)dBd(t> (18) 
Yftf, e)Yp\cl>, 6) S m(<t>)ded<t> = 5 u ,5 kk ,. 



o J o 



J 



Applying Theorem [4] to the system {Q^} , whose sampling matrix Q^(4>j,9j) = yj is equivalent 
to the preconditioned spherical harmonic matrix (sin0j) 1 / 2 l^ s (0y, 9j) = Uj Theorem 1 follows 
from the recovery results for restricted isometry systems in Theorem |3j 
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